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QCD sum rules involving mixed inverse moment integration kernels are used in order to determine 
the running charm-quark mass in the MS scheme. Both the high and the low energy expansion 
of the vector current correlator are involved in this determination. The optimal integration kernel 
turns out to be of the form p{s) = 1 — (so/s)^, where so is the onset of perturbative QCD. This 
kernel enhances the contribution of the well known narrow resonances, and reduces the impact 
of the data in the range s ~ 20 — 25 GeV^. This feature leads to a substantial reduction in the 
sensitivity of the results to changes in so, as well as to a much reduced impact of the experimental 
uncertainties in the higher resonance region. The value obtained for the charm-quark mass in the 
MS scheme at a scale of 3 GeV is mc(3 GeV) = 987 ± 9 MeV, where the error includes all sources 
of uncertainties added in quadrature. 



PACS numbers: 12.38.Lg, 11.55.Hx, 12.38.Bx, 14.65.Dv 

Progress on the theoretical |lj-|13|. as well as on the 
experimental information [l3j-[l6j| on the vector current 
correlator has allowed for a considerable improvement 
on the accuracy of QCD sum rule determinations of 
the charm-quark mass [13- [H. The analysis of ^ 
is based on inverse (Hilbert) moment QCD sum rules, 
requiring QCD knowledge of the vector correlator in 
the low energy, as well as in the high energy region. In 
[l8| an alternative approach was used which involves 
only QCD information at short distances, together 
with (a) a simple integration kernel p{s) = 1 — s/so 
(local constraint), and (b) Legendre-type polynomial 
kernels (global constraint). In this paper we describe an 
improved analysis based on the use of direct as well as 
inverse moment kernels of the form p{s) = 1 — (sq/s)"', 
with n > 1. These kernels enhance considerably 
the impact of the well known narrow resonances, as 
compared with e.g. a simple kernel p{s) — 1/s^ , or 
p{s) = 1 — s/so- They also provide a welcome stronger 
suppression of the contribution of data in the range 
s ~ 20 — 25 GeV^. In comparison with simple inverse 
moments without pinching, this means that results 
are less sensitive to assumptions about the onset of 
perturbative QCD (PQCD), as well as to the treatment 
of the higher resonance data. For instance, changes in 
So in the range Sq — 15 — 23 GeV'^ lead to a variation in 
TOc(3 GeV) of only 4 MeV (for n = 2) as opposed to a 
variation of 14 MeV for p{s) = 1/s^, as used in \17^ 



We consider the vector current correlator 

n^.(9') = i I d^a;e'«^(0|T(F^(x) K(0))|0) 



(9p (ii^ - q'^dt^iy) n((?^) , 
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where V^{x) = c(x)7pc(x). From the residue theorem in 
the complex s-plane {—q^ = = s) it follows 
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p{s)n{s)ds 



+ Res[n(s)p(s),s = 0] , (2) 
where p{s) is an integration kernel, and 
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with Rc{s) the standard R-ratio for charm production. 
The PQCD expansion of n(s) at short distances can be 
written as 

n(.)|PQcz. = e^E(^)"n(")(.), (4) 
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where ec = 2/3 is the charm-quark electric charge, and 



=0 



n. 
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and TOc = ^c(m) is the running charm-quark mass in the 
AfS'-scheme. Up to order O [al{fhl/s)^] the function 
n(s)pQcr> has been calculated in and exact results 
for IIq'^^ and H^'' have been found in 0. The function 
n'^^ is known exactly up to a constant At five-loop 
order ©(af) the full logarithmic terms for IIq'*' may be 

found in [5] , and for 11^^^ in . Since there is incomplete 
knowledge at this loop-order we shall use the available in- 
formation as a measure of the truncation error in PQCD. 
The low energy expansion of the vector correlator around 
s — can be written as 



3ei_ 
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n>0 



(6) 



where z = s/(4m^). The coefficients C„ can be expanded 
in powers of as(^) 



/^(30) , ^(31), , ^(32),2 



(7) 



where Im = ln(TO^(^)/^^). Up to three loop level the 
coefficients up to n = 30 of €„ are known J8[- p. At 
four- loop level we have Co and Ci from [8] , [10| , C2 from 
0, and 6*3 from [Tl'|. We wiU choose p{s) so that no 
coefficients C4 and above contribute to the residue at 
s = 0, Res[n(s)p(s),s = 0]. 

Apart from the quark mass, the fundamental QCD pa- 
rameters are the running strong coupling as(M^), and 
the gluon condensate. For the strong coupling we use 
the world av erag e from [23|, which agrees with lattice 
QCD results [2lf, as(M|) = 0.1184 ± 0.0007. However, 
we will consider other values when comparing results for 
TOc with other analyses. In the non-perturbative sec- 
tor the leading power correction in the Operator Prod- 
uct Expansion (OPE) involves the gluon condensate, i.e. 
((as/7r)G^) whose value has been extracted [13] from 
the ALEPH data on r-decays. While the gluon conden- 
sate is renormalization group invariant, its determina- 
tion from QCD sum rules involves a difference between 
integrals of PQCD and integrated experimental data. 
This leads to an unavoidable dependence of the gluon 




FIG. 1: Experimental data for the total R{s) ratio together 
with the optimal integration kernel, Eq.(8), with n = 2 (dash 
curve), and p{s) = (solid curve) normalized to coincide 
with the former at the position of the ip{2S) peak. 



condensate on the value of as used in the PQCD ex- 
pression of the correlator. Extrapolating the results of 
|23 | to include current values of [1^, 23 1, leads to 
{{as/-K)G^) = (0.01 ± 0.01) GeV^ This large uncer- 



tainty in the value of the gluon condensate has only a 
very small impact on our results for rn-c- 
Turning to the experimental data, we follow closely the 
analysis of [Tt} . For the first two narrow resonances we 
use the latest data from the Particle Data Group p^ ]. 



M, 

M, 



3.096916(11) GeV, P 
= 3.68609(4) GeV, P 



j/V^e+e- = 5.55(14) keV, 
V>(2s)^e+e- = 2.35(4) keV. 
These two narrow resonances are followed by the open 
charm region where the contribution from the light quark 
sector Ruds needs to be subtracted from the total R-ratio 
Rtot- We perform this subtraction as in [2^. In the re- 
gion 3.97 GeV < < 4.26 GeV we only use CLEO data 
|16| as they are the most precise. In connection with the 
three data sets from BES 13|-[l^, we assume that the 



systematic uncertainties are not fully independent and 
add them linearly, rather than in quadrature. However, 
we treat these data as independent from the CLEO data 
set fl^, and thus add errors in quadrature. There is no 
data in the region s = 25 — 49 GeV^, and beyond there 
is CLEO data up to s ~ 110 GeV^ The latter data is 
fully compatible with PQCD. 

We discuss next the integration kernels p(s) in Eq.(2), 
which we choose as 



(8) 



with n > 1. As discussed in 17|, inverse moments 
p(s) — 1/s" should not involve too large values of n. 
In fact, the convergence of PQCD deteriorates with 
increasing n, the gluon condensate contribution increases 
sharply for n > 2, and the uncertainties in as and the 
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Kernel 












1129 


1021 


998 


995 


1 - {so/sf 


1146 


1019 


991 


987 



TABLE I: Results for the charm-quark mass at different orders in 
PQCD, and for two integration kernels. The results for p(s) = 
are obtained using slightly different values of the QCD parameters, 
and a different integration procedure as in [TtII . 







Uneertainties (in MeV) 




Kernel me(3GeV) 


EXP 


Aoa A/J NP Sq 


Total 


s^^ 995 


9 


3 1 1 14 


17 


l-{so/sf 987 


7 


4 114 


9 



TABLE IL The various uncertainties due to the data (EXP), the 
value of as (Aos), changes of ±35% in the renormalization scale 
around fi = 3 GeV (A/i), the value of the gluon condensate (NP), 
and due to variations in so (so)- 



renormalization scale /i have a greater impact on the 
total error of the charm-quark mass. On the other hand, 
direct kernels of the form p{s) = s", with n > 1, pose 
problems. Indeed, the high energy expansion of the 
vector correlator is incompletely known at ©[a^], so that 
the greater the value of n the greater is the contribution 
of the higher order mass corrections at this order in 
PQCD. Already with n = 1 one would need a Pade 
approximation for the term ©[a^m^]. Hence, in order 
to avoid any approximation up to this order in PQCD 
we restrict ourselves to direct moments with n = 0, and 
include inverse powers in an attempt to enhance the 
contribution of the well known narrow resonances, J/V' 
and ■0(25'), and at the same time suppressing the broad 
resonance region. We found that Eq.(8) with n = 2 
is the optimal kernel as explained next. In Fig.l we 
show the experimental data for the ratio R{s) together 
with the kernel Eq.(8) with n = 2 and sq ^ 23 GeV^, 
and the simple kernel p{s) — 1/ normalized such that 
both kernels coincide at the peak of the second narrow 
resonance -0(25'), i.e. s ~ 13.6 GeV^. One can easily 
appreciate that in comparison with the latter, the former 
kernel leads to a welcome enhancement of the weight of 
the J/'0, as well as to a strong suppression of the broad 
resonance region, and particularly the region near the 
onset of the continuum. Quantitatively, the ratio of the 
area under the hadronic spectral function weighted with 
p(s), in the narrow resonance region, and in the 

broad resonance region and beyond, T^^, is I„r/2^br = 3.6 
forp(s) — 1/s", and I„r/2^br = 7.7 for Eq.(8) with n = 2. 
Other values of n lead to slightly less enhancement. In 
addition, the kernel Eq.(8) with n — 2 leads to final 
results for rric which are fairly insensitive to the choice 
of So- For instance, in the range Sq — 15.0 — 23.0 GeV^, 
TOc(3 GeV) changes by about 4.0 MeV, while using the 
kernel p{s) = it changes by 14.0 MeV. 
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FIG. 2: Results for mc(3 GeV) as a function of so for the 
kernel p(s) — The variation of rnc(3 GeV) in this range 

is up to 14 MeV. 
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FIG. 3: Results for mc(3 GeV) as a function of so for the 
kernel Eq.(8) with n = 2. The variation of mc(3 GeV) in this 
range is up to 4 MeV. 



Proceeding to our determination we list in Table 1 the 
results for mc(3 GeV) at different orders in perturbation 
theory, and using two integration kernels. The results 
for the kernel p{s) — differ from [T^l as we use now 
a slightly different value of the strong coupling, of the 
gluon condensate, and of the ip{2S) parameters as given 
above, and we include the CLEO data [l^. The vari- 
ous errors associated to the final value of mc(3 GeV) are 
given in Table 2. Results from Fixed Order Perturbation 
Theory (FOPT) are essentially the same as using Con- 
tour Improved Perturbation Theory (CIPT) to integrate 
around the circle of radius sq in the complex s-plane. In 
Fig. 2 we show the results for TOc(3 GeV) as a function 
of -So for the kernel p{s) = and in Fig. 3 for the 

kernel Eq.(8) with n = 2, the latter exhibiting improved 
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stability. 

The convergence pattern in as of the PQCD integral as 
a function of rric can be studied by computing 

^{so) = --^ / pis)U{s)ds+Res[U{s)p{s),s^ 0] , 
2^^* Jc(|.ol) 

(9) 

with these integrals being functions of both rric and 
a^. Using 771^(3 GeV) = 987 MeV, and Eq.(8) 
with n = 2 we find reasonable convergence, i.e. 
T(°) = 91.4 GeV^ 1(1) = 62.0 GeV^ I^^) ^ 57 q GeV^ 
and 1(3) = 56.3 GeV^ where the upper index in X^^^ 
indicates the power of a^. 



Our final result using the optimal kernel, Eq.(8), with 



n = 2 is 

mc(3 GeV) = 987 ± 9 MeV , (10) 

in good agreement within errors with the result from in- 
verse mome nt Q CD sum rules [l7j , other recent determi- 
nations [ll|-[i3, m, nil, as well as lattice QCD dJ. 
Translated into a scale invariant mass, the above result 
gives ffic{fhc) = 1278 ± 9 MeV for the value used here 
for the strong coupling. 
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